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FLOW  FIELD  ON  THE  LEE 
SIDE  OF  A  DELTA  WING 


I.  Introduction 


This  study  is  concerned  with  the  flow  field  on  the 
lee  side  of  a  flat  plate  delta  wing  with  supersonic  leading 
edges.  As  shown  in  Fig.l,  the  wing  is  placed  at  an  angle 
of  attack  a  and  at  zero  yaw. 


The  vertex  of  the  wing  is  placed  at  the  origin  of  a  rec¬ 
tangular  Cartesian  coordinate  system.  The  flow  medium  is 
assumed  to  be  an  inviscid,  non-heatconducting,  nonreacting, 

Preceding  page  blank 
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Ideal  gas.  The  free  stress  Is  assumed  to  be  uniform  end 
steady . 

The  general  flow  pattern  on  the  expansion  side  of  the 
wing  Is  symmetrical  with  respect  to  the  XZ  plane  and  ? 
sketch  of  the  flow  pattern  on  the  lee  side  Is  given  in 
Fig.  2. 


Figure  2.  Plow  Pattern  on  Surface 

As  the  free  stream  flows  over  the  leading  edge,  the 
component  of  velocity  parallel  to  the  leading  edge  is  not 
altered  but  the  component  normal  to  the  leading  edge  must 
expand  to  become  parallel  to  the  surface  of  the  wing. 

This  results  in  the  total  velocity  vector  on  the  surface  of 
the  wing  immediately  downstream  o^  the  leading  edge  having 
a  component  normal  to  the  root  chord.  Hov:ever,  because  of 
the  symmetry,  the  total  velocity  vector  at  the  root  chord 


must  he  oarellel  to  the  root  chord 


mere  fore,  as  indicated 
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in  Pig.  2,  the  total  velocity  vector  on  the  surface  of  the 
King  nust  change  direction  and  this  condition  is  brought 
about  by  the  formation  of  an  internal  shock  wave  which  is 
attached  to  the  upper  surface. 

The  flow  field  for  a  flat  plate  delta  wing  with  super¬ 
sonic  leading  edges  is  conical  as  conceived  by  Busemann 
(Hef  5)  and,  as  such,  the  flow  field  in  any  plane  normal  to 
the  root  chord  has  a  similar  solution.  Thus,  the  entire 


lee  side  flow  field  may  be  represented  by  any  arbitrarily 
selected  plane  which  is  normal  to  the  root  chord.  Tnls 
plane  Is  called  the  "crossflow”  plane  and  a  qualitative 
sketch  of  the  right  half  of  the  plane  is  shown  in  Pig.  3* 
The  crossflow  plane  may  be  separated  into  six  regions 
with  distinct  characteristics  as  indicated  in  Fig.  3.  The 
goverr.inr  equations  in  regions  I,  1?.,  and  III  are 
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hyperbolic.  The  governing  equations  in  regions  F/,  V,  ai*d 
VI  are  elliptic.  Line  3D  represents  a  crossflow  scnic  line 
along  which  the  governing  equations  are  parabolic.  Line  DSG 
is  an  internal  shock  wave.  The  segment  GE  is  planar  but  the 
segment  DS  is  curved.  Line  A3  represents  the  first  ray  of  a 
Prandtl-Meyer  expansion  fan.  Line  3FE  is  the  reflection  of 
this  expansion  ray  froa  tne  surface  of  the  cone  of  excita¬ 
tion  BC  which  emanates  froa  the  vertex  of  the  wing.  Line  A? 
represents  the  last  ray  of  the  Prandtl-Keyer  expansion  fan. 
Lines  CD  and  OS  are  crossflow  streamlines  separating  the 
rotational  flow  in  region  V  froa  the  irrotational  flow  in 
regions  IV  and  71.  line  OGA  is  the  wing  surface  and  line  OC 
represents  the  plane  of  symmetry. 

Among  the  analytical  solutions  to  this  flow  field  are 
those  presented  by  Haslen  (Bef  13) »  Babayev  (Bef  1).  3eesan 
and  Powers  (Ref  3)»  and  Kutler  and  Losax  (Bef  li).  Fowell 
(Ref  8)  presented  an  analytical  solution  along  with  experi¬ 
mental  data  for  the  flow  field.  Bannlnk,  Kebbllng,  and 
Reyn  (Ref  2)  conducted  an  experimental  Investigation  of  the 
flow  field  as  did  Cross  and  Hankey  (Ref  7).  A  compilation 
of  papers  treating  the  mathematical  theory  of  conical  flovs 
is  presented  in  the  work  of  Bulakh  and  Reyn  (Ref  4). 

In  Kaslen’s  treatment  of  the  problem,  it  was  assured 
there  were  no  internal  shock  waves  in  the  flow  f'iela  and  a 
solution  was  obtained  by  using  Prandtl-Meyer  relations  in 
regions  I  and  II  and  a  relaxation  method  in  the  remainder  of 
the  flow  field.  Fowell  experimentally  determined  the 
existence  of  an  internal  shock  wave  and  then  assured  the 
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shoes  wave  was  planar  and  used  an  approximate  method  to 
locate  the  shock  me  and  solve  for  the  flow  adjacent  to 
the  win*  surface,  The  solution  presented  by  Babayev  also 
included  an  internal  planar  shock  wave  but  the  flow  field 
in  regions  III,  17,  V,  and  VI  was  determined  from  irrota- 
tional  potential  flow  theory  using  the  method  of  successive 
approximations.  Seeaann  and  Powers  assumed  the  existence  of 
an  internal  shock  wave  and  attempted  to  solve  the  problem  by 
using  a  modified  method  of  characteristics  scheme.  However, 
as  a  result  of  difficulties  with  the  Internal  shock  cal¬ 
culations,  only  shock-free  results  were  reported.  In  the 
experimental  work  of  Bannihk,  Hebbeliug,  and  2eyn,  the 
existence  of  an  internal  shock  wave  was  reconfirmed  but  the 
shock  wave  was  found  to  be  curved. 

In  all  of  the  analytical  studies  cited  above  it  is 
assumed  that  flow  in  region  III  is  dependent  on  the  flow  In 
regions  IV,  V,  and  VI.  The  approach  taken  in  this  study 
Is  to  solve  the  flow  field  in  region  III  by  reducing  the 
lnviscid  flow  governing  equations  using  conical  flow  con¬ 
ditions  and  then  applying  the  method  of  characteristics  to 
the  resulting  hyperbolic  equations.  The  results  of  this 
study  may  then  be  applied  to  the  solution  of  the  flow  field 
in  the  elliptic  region. 
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II.  The  Development  of  Equations 

The  crossflow  plane  say  fee  taken  as  any  YZ  plane  that 
intersects  the  positive  X-axis.  This  plane  is  noniisen- 
sionalized  by  using  the  dimensionless  coordinates  and  , 

Z 

where  £  =  -  (1) 

A 

v 

and  ’  -  -  (2) 

X 

The  governing  equations  are  transf orsed  to  the  crossflow 
plane  by  seans  of  these  dimensionless  coordinates. 

Line  BC  -  Cone  of  Excitation 

The  cone  of  excitation  which  Is  generated  by  the  apex 
of  the  delta  wing  Is  given  by, 

r~2  2 

yY  +  (-X  sin  a  +  Z  cos  c) 

-  (X  cos  a  +  Z  sin  a)  tan  u=  0  (3) 

In  terms  of  i  and  n  this  becomes 

j— - ^ 

y n  +  (icos  a  -  sin  a)  -  (cos  a  +  i  sin  a)  tan  0  (4) 

or,  since,  tan  ^  =  (K  -  1)  ' 
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Cl - ?-)  (l  i*  *  v2)  =  (cos  a  *  £  sin  c  )2  (5) 


3eglon  ^  -  Prandt I -Sever  Szuar.slon  Pan 

The  development  of  the  equations  in  this  section 
follows  closely  that  of  rabayev. 

Fron  the  classical  Prandtl-Keyer  relationship,  as 
given  in  liepsann  and  Hosh>o  (3ef  12:93-100),  the  local  flow 
deflection  i  in  a  plane  norsal  to  the  leading  edge  is 
given  by 


or,  by  letting 


-1 


e  =  4' -  tan  *  v 

T  Y-l  \ 


fFI 


v+5 


n 


(7) 


and 


v  =  tan 


-1  rr 


r»  —  « 
!  n 


(3) 


equation  (6)  reduces  to 


6=0-  tan 


_  1 


-r  . 


(9) 


w*.  cr 
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The  direr.slonless  velocity  components  In  region  I  are 
c! men  as  functions  of  6  by 


sa  -  fe  cos  (e^-6)  cos  x  +  cos  £  sin  x 


v  =  -o  cos  (c--5)  sin  x  +  cos  £  cos x 


(10) 


(II) 


xnere 


w  =  b  sin  (.>~6) 


aj  =  tan 


j  /  tan  a  \ 

l  cos  x  I 


(12) 


U3) 


£  =  cos~  *  (cos  a  sin  x  ) 


(1*) 


fe2  =  -L  (1-  cos2^  ) 

Y+l  4  Y+1 


r  2  2  pi 

x  + -  sm  * ' -  a 

Y-l  iV1 


(15) 


Equations  (10),  (11),  ana  (12)  are  developed  In  detail  In 
Appendix  A. 

In  the  plane  normal  to  the  leading  edge  the  angle 
between  a  Prendtl-Keyer  expansion  ray  and  the  surface  of  the 
wing  is  given  by 


o  =  e  -  +  —  -  0  +  a1 

oc  00  2  1 


0  6) 


By  noting  that 


tan  c  = 


cos  >  -  i  sm  x 
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the  slope  of  the  expansion  ray  In  the  £ ,  t  plane  say  he 
obtained  by  sabscitutin?  equation  (17)  Into  equation  (l6) 
and  differentiating.  This  results  in 


exp  ray 


n 

-  sin  x  tan  (  ©w-  +  — 


a  _  S  ) 
1 


(131 


setting  £  =  0  In  equation  (17)  gives  the  * -value  of  the 
Intersection  of  all  expansion  rays  with  the  !t-axis  as 


*  =  cot  x 

sax 


(19) 


Sanations  (iu),  (ii),  ( 12),  (18),  am  (19)  completely 
determine  the  dimensionless  velocity  distribution  through¬ 
out  region  I.  Other  flow  properties  may  be  determined  from 
the  Isentropic  relationships. 


Beglon  II  -  Fully  Expanded  Area 

fhe  flow  properties  In  region  II  are  determined  com¬ 
pletely  by  the  values  along  the  last  expansion  ray.  This 
is  a  result  of  the  fact  that  the  flow  properties  are  con¬ 
stant  throughout  the  region. 

Setting  6  =  in  equation  (9)  results  In  the 
expression  for  the  value  of  6  along  the  last  expansion  ray 
as 


e 


cr 


+  6co  -  -CO  + 


pFT 
tan  ..  - 

>  vn 


(20) 


9 


ni*r 


Point  3  -  Intersection  of  Cone  of  Excitation  ana  First 
Expansion  3ay 

The  cone  of  excitation  Intersects  thv  first  expansion 
ray  at  the  point  Mere  the  slopes  of  these  lines  are  equal. 
From  equation  (13)  the  slope  of  the  first  expansion  ray  Is 


% 

d  t  \  sin  x 

d*  /exp  ray  1  tan  (  vk-  £>-  ) 


(21) 


and  differentiation  of  equation  (5)  results  in 


d’ s  /cone 


(i - _), 


(1 — y~)  i  -  sin  c  (cos  c  + 
iC 


c  sin  e.  ) 


(22) 


equating  equations  (2i)  and  (22)  results  in  the  i  inter¬ 
section  value  i .  .  of 

X  ilv 


i  tan  (v^-a.) 

(1  -  -«-)  - ; - =~  +  sin  c.  cos  a 


sin  x  tan  x 
f  /  . 


(1  - 


1  /  /  ./  _  %  I  C-  % 

^  \ 1+  nrr  “  i 
_  _ 


-  sin  a 


The  intersection  value  ^nt  nay  be  obtained  fron  equation 

(5). 
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Seglon  III  -  Pseu-io- £111  otic  Area 

The  deTelopcerst  of  the  eqaatlcrs  for  r-.*&ion  III 
closely  follows  that  of  Chains  and  i'sgaer  (S;f  6)  vho 
present  a  detailed  derivation  of  tee  governing  equations  for 
raon-isentroplc  flex. 

fhe  continuity,  aos-entus,  and  energy  equations  can  be 
written  In  vector  fora  as 


V  •  (  p  U)  =  0 

(24) 

ff  I  U  =  x  vs 

(25) 

F  •  vs  =  o 

(26) 

Using  the  conical  flow  conditions,  the  equations  aay  be 
rewritten  for  the  i ,  ”  plant*  as 


d'J 

2A  (rj  -  + 


du  gv  d» 

d*  bn 


-  (nu  -  v) 


(u2  +  v"  +  w2  ) 
bn 


-  (iu 

-  w) 

b 

(u2  +  V2  + 

w2  ) 

=  0 

bk 

av 

du 

aw 

aw 

au 

v  (  *>  - 

+  $  — 

+  - ) 

+  W 

(  n  -  +  t 

a* 

6sc 

bn 

br, 

bk 

bk 

as 

—  - 

-  A  (  n 

—  + 

(  — ) 

6*1 

bk 

(Zb  1 
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u 


B 


a*  a*  a*  a55 

*  _  +  - — )  +  v  (- —  — )  =  -  a  — - 

a$  a?  ^  a$  a? 

a*  a^  .  a^  a*  %  a« 

£  -  -r  - J  +  V  ( -  -  - )  =  -  A  - 

a*  a*  a*  a*  a$ 


as  as 

(n U-Y )  -  +  (fu-w)  -  =  0 

a*  a? 


(29) 


(30) 


(3i) 


where  A  is  given  by 


i  Y-l  2  2  2. 

A  =  +  -  (i-u  -v  -if  ) 

►£,  2 


(32) 


Equation  (23)  is  a  linear  combination  of  equations 
(29).  (3°).  and  (31)  and  it  may  then  be  discarded.  By 
manipulating  equations  (27)*  (29)*  (30).  ana  (31)  and  using 


the  definition  of  a  total  derivative  the  following  set  of 
equations  may  be  written  in  matrix  form  as 


a 

1 

a2 

a3 

fil 

£> 

a? 

G 

1 

0 

0 

-1 

av 

<£> 

Q 

d* 

dv 

0 

0 

(r’ 

a? 

dv 

0 

0 

d( 

dv 

aw 

(-) 

a»j 

dv; 

(33) 
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xaere 


a.  =  Aiv  -  {iu-v)  (£u-w) 


(3*0 


a?  =  A  (l+n  )  -  (iju-v)4 


(35) 


a3  =  A  (l+{2)  -  (iu-w)2 


(36) 


_  A  ds  ♦  u  (du  +  »»dv  +  £dw) 
(»ju-v)  d£  -  (|u-w)  dn 


(37) 


A  ((u  -  >jw)(du  +  ndv  +  £dw) 
(iu-v)  d£  -  (£u-w)  dv 


(38) 


Setting  the  coefficient  matrix  of  equation  (33)  equal 
to  zero  gives  the  slope  of  the  characteristics  as 


al±  ^-a2a3 


d1?  /  ±c 


(39) 


The  compatibility  equations  for  these  characteristics 
are  obtained  by  setting  the  determinant 


al  a2  a3  G 


0  0  Q 
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Substituting  equations  (37)  and  (38)  into  equation  (40) 
results  in 


where 


-  b^du  +  bgdv  +  b^dw  +  b^ds  =  0 


b^  =  -A  (£  v-nw) 


b2  =  -An(^v-j?w) 


-[va2<Jo±]“ 

-  [al*a2(~)o±]’,U 


£  ({u-l 


+  a2  I  (£u-w)  -  ( 7j  u-v )  ( 


a* 

d^c±  . 


b^  =  -A^(^v-nw)  - 


d* 


t  u 


r 

*3  ($U-w)  -  (>?U-v)  (— •) 

J  aij 

L 


C± 


d£ 

(-) 

d*? 


c± 


r  d* 

D/,  =  -A  ai  ~ap  ( —  '  + 

L  *  dn  c  -  _ 


(41) 

(42) 


(43) 


(44) 


(45) 


Since  the  flow  properties  in  the  free  stream  are  taken  as 
uniform  and  there  are  no  rotation  producing  effects  in 
regions  I,  II,  and  III,  It  is  reasonable  to  assume  the  flow 
in  the  region  III  is  isentropic.  Equation  (4l)  then  reduces 
to 

b^du  +  b2dv  +  b-jdw  =  0 
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For  any  streamline 


U  •  n  =  0  (4?) 

where  n  is  a  unit  vector  normal  to  the  streamline.  In  the 
crossflow  plane  this  results  in  the  slope  of  the  streamline 
being  given  as 


£u-w 

7JU-V 


(48) 


A  compatability  equation  for  the  streamline  may  be  obtained 
by  multiplying  equation  (29)  by  (nu-v)  and  equation  (30)  by 
(£u-w)  and  adding  them.  This  results  in 


du  dv  dw 


d» 


d^ 


dw 


(r,u-v){  —  +1J—  +£--  |+  (*u-w)  1=0 


d>?  d*/  d»j 


at 


(49) 


Substituting  equation  (48)  and  using  the  definition  of  a 
total  derivative  reduces  equation  (49)  to 


du  +  *j  dv  +  £  dw  =  0 


(50) 


which  is  valid  only  along  a  streamline. 

Line  BD  -  Crossflow  Sonic  Line 

Since  the  governing  equations  are  parabolic  along  the 
sonic  line,  the  term  under  the  radical  in  equation  (39)  must 
be  zero  at  a  sonic  point.  This  results  in  both  character¬ 
istics  having  the  same  slope  at  a  sonic  point  and  the  two 
compatability  equations  represented  by  equation  (39) 
becomes 
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(50 


IB*  l®»s  of  ©a*  of  the  ©gsrp&fcl  ©1 11  ty  reiatl>E3higKS  Is 
caasgessated  for  by  the  fast  that  the  ercssflcw  Sach  nsarber 
E  acsst  be  osnity  at  the  sonic  Use,  IB*  crossflow  Each 

<F 

masher  Is  given  by 


Boundary  Conditions 

IBs©  equations  obtained  by  writing  equation  (*rl)  for 
both  characteristics  am  equation  (5C)  for  the  streaaliEe 
provide  a  set  of  three  ordinary  differential  equations  for 
u,  ?,  and  w  which,  with  sui  table  boundary  conditions,  ear 
be  used  to  solve  for  the  dirensionless  velocity  distribution 
in  region  III.  A  set  of  suitable  boundary  conditions  for 
this  area  consists  of  the  dirensionless  velocity  distribu¬ 
tion  along  line  BF2  and  the  condition  that  the  crossflow 
velocity  is  sonic  along  line  BD^ 
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Ill.  'Pistracta*  1  cn  35c® 


Seller;  I  —  PranefE— Serer  -•gr^-r^lsr.  Fas 

©sirs  laevt® sBs  methrt&i  of  taraxeirts,  «iosatt©E  (2®)  Is 
sclrei  for  ?cr-  'Ihe  slc-ce  of  IrilTldtsal  expansion  rajs  Is 
then  obtained  frem  es-atior  (l3)  for  eo^al  Increments  of  5 
Is  the  range  from  to  5  amd  nsec  art. 42s  equation  (19)  to 
establish  the  loess  of  the  expansion  rajs,  ibe  Telocity 
ccocporerts  alor.-r  each  raj  are  found  by  solxlng  equations 
(9)  arid  (15)*  using  the  respective  values  of  5,  and  sub¬ 
stituting  the  resalts  I sto  equations  (10),  (11 ),  and  (12). 


fine  oFE  -  Eefiecilon  of  First  scarsioa  Say  Frege  the  Sonic 
l  i  y-^=- 

rhe  locas  of  the  reflection  of  the  first  expansion  raj 
is  found  by  an  Iteration  technique  in  which  the  initial 


assumption  is  that  the  reflection  has  a  vertica  l  slope  be¬ 
tween  the  points  B  arsd  A^  which  are  on  the  successive  expan¬ 
sion  rays  i  and  i+i  as  shown  in  Fig.  4.  The  position  of 
point  B  is  known  since  It  is  the  intersection  of  ray  1  and 
the  first  reflection.  Point  is  the  intersection  of  ex¬ 
pansion  ray  i+i  ana  the  first  reflection  ana  is  detennineH 
by  equations  (23),  (5),  (l8),  ana  (15). 

The  reflection  of  the  first  expansion  ray  is  a  charac¬ 


teristic  for  region  III  and 


* ts  slope 


ac 


is 
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©btatl oei  by  substituting  5.  ,  *  ,  and  ftfee  Telocity 

~I  *1 


COTraoneEts  for  ray  {i*lJ  Ieco  equatlos  (39).  ISslcg  this 

valsse  of  slope  a  straight  Use  Is  passed  through  point  jt  _ 

1 

Urf?  I  e  tors  set*  m*  c  of  this  line  rsy  i  Is  f ©mass , 


distance  between  C  and  5  Is  then  determined  and.  If  It  Is 
larger  than  specified  Units,  tine  procedure  Is  restarted  by 
assuming  the  value  of  d{\  at  for  the  slope  of  the 

a7?  /c+ 

first  reflection  at  point  5  and  determining  point  as  the 
new  intersection  of  ray  1+1  and  the  first  reflection. 

Lire  BD  -  Crossflow  Sonic  line 


A  sketch  indicating  the  geometry  used  in  solving  for 
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<f’&rre  A3  Is  ths-  ccat  ®f  azcitamlcm.;  liases  5C  and  Dc.  ere 
ezoisaESl®«  rays;  its*  53  Is  a  segment  ©f  the  reflect  i.(D*2S  (&!T 
the  first  esrarsicas  ray;  lines  31  and  31*  are  first  family 
characteristics ;  arc  lines  5F  33  are  streamlines. 

Alt hoiBgSs  the  sectnetric  sletch  in  Pig.  5  is  related  to 


Figure  Locating  Sonic  Point 


solution  of  the  first  sonic  point,  the  procedure  provided 
below  is  valid  for  all  sonic  points. 

The  streamline  slope  based  on  the  properties  at  B  is 
computed  using  equation  (43)  and  a  line  Kith  this  slope  is 
passed  through  5.  The  slope  of  the  first  family  character 


istic  based  on  the  "low  conditions  at  D 


is  computed 


from  equation  (39)  and  a  line  with  this  slope  is  passed 
through  D  to  intersect  the  streamline  at  P. 
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*i\ 

Hhe  <3 if  —  I  Is  nabs  titrated  la  tSss  left  side  ©f 

cquiattcm  C5il-  Th*  re-smiting  eqtsatlom  Is  used  with  eqma- 
fci©a»  (46)  sacf  {5®?  t©  s©Xre  for  fee  velocity  ©©mpememts 
at  ?. 

ffee  slope  of  the  streamline  at  ?  Is  then  ©©mpsted  from 
e©T!£S±I©m  (43)  ami  averaged  Vi  fe  the  slope  of  the  streamline 
at  E.  A  streamline  arife  this  average  slope  Is  passed 
through  5  to  Intersect  characteristic  3a  at  G.  «ith  fee 
value  of  fee  slope  of  Use  DS  substituted  into  equatioE 
(51),  equations  (to),  (5®)*  and  (51)  are  solved  for  the 
Telocity  components  at  G.  the  slope  of  fee  streamline  at  G 
is  then  computed  from  equation  (43)  and  averaged  xith  the 
slope  of  line  53.  A  streamline  with  this  averaged  slope  Is 
passed  through  B  tc  intersect  characteristic  BI  at  H.  with 
the  value  cf  the  slope  of  line  DE  substituted  into  equation 
(51 )#  equations  (46),  (50),  and  (5l)  are  solved  for  the 
velocity  components  at  H.  These  velocity  components  and  the 
coordinates  of  H  are  then  substituted  into  equation  (52)  to 
determine  the  crossflow  Kach  number  at  H. 

If  the  crossflow  Mach  number  is  not  sufficiently  close 
to  unity,  a  r.ew  point  H‘,  which  is  located  at  a  small 
distance  from  H  and  on  the  streamline  BH,  is  selected. 

Using  the  slope  of  the  line  DH* ,  the  velocity  components  at 
H1  are  determined  and  the  crossflow  Kach  number  computed. 
This  procedure  is  repeated  until  a  position  on  the  line  BH 
is  found  where  the  crossflow  Koch  number  is  x.'ithin  ±0.001 
of  unity. 
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if  (ter  £  scseLc  ?©*r(t  ha©  beer  fssai.  the  float  ©©adltlccs 
al©2jg  the  secozti  family  charactcerlrtl©  ssssZjss  Shr®suxh  the 
r©nl©  ?<2>lrtt  are  iettermlasd  by  (tee  ©rccedur©  described  in 
the  fcllcwlng  section.  Saawe^er,  if  the  positlcs  ©f  any 
point  al©*3g  this  characteristic  ©an  net  be  determined  within 
the  specified  limits,  a  cxr*  p©Int  2s  on  the  line  23  is 
selected  £s  the  point  throsigh  which  t©  pass  the  initial 
streamline  ard  the  entire  scnic  point  pro©*© tire  is  repeated. 

region  HI  —  Field  golnts 

£  geometric  sketch  relate©  t©  the  solution  ©f  the  flow 
field  in  region  III  Is  given  is  Fig.  6.  Although  the  sketch 
ami  the  procedure  outlined  below  refer  to  a  particular  point 
In  the  region,  the  procedure  is  completely  general,  fhe 
lines  lab  led  by  I*s  and  II Ss  are  first  and  second  family 
characteristics  respect! rely.  Uses  ED  and  2*D*  are 
streamlines. 

Using  the  known  conditions  at  3  the  slope  of  the 
streamline  is  computed  from  equation  (43)  and  a  line  with 
that  slope  Is  passed  through  B.  Using  the  conditions  at  C 
the  slope  of  the  first  family  characteristic  is  computed 
and  a  line  with  this  slope  is  passed  through  C  to  intersect 
the  streamline  at  D.  In  this  case,  as  opposed  to  the  sonic 
line  ease,  equations  (*Vi)  represent  two  characteristic 
compatability  relations  and  these  relations  are  used  with 
equation  (5C)  to  solve  for  the  velocity  components  at  D. 

The  conditions  at  D  are  then  used  to  compute  the  slope 


nanrer  a."  b. 


The  distance  S’  to  A  is  detcrnir.ee  and,  if  i 
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is  rao'S  vitals,  limits,  ?/b*»  eriire  srsc-tss  *  s  re?«att«rf  acltfe 
s  sear  jeinS  3*. 


2~< 


Ix.  SS-sceas*  os.  of  Sa-rslts 

Saner  los!  solr.tjorcs  were  oSrialiaei  far  t&s  foiloaflsr 
sets  of  Isoat  oata  vlth  1ZG  exz&ns Ion  rays  la  resist  Is 

1)  =  2.56,  a  =  1^.2®,  x  =  £5.3®  ; 

2 )  ft.  =  %.©,  a  =  5-®°.  *=  6©.C°  ; 

3)  ^  =  3-3,  a  =  i2.©°,  x  =  £5.©°  . 


in  ‘.'■.e 


Figure  7.  Characteristic  Network  for  fi>«*2.96 
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ares,  of  the  Intersection  ef  the  X&ch  cone  arj£  che  first 
expansion  raj.  Bcxerer,  the  amount  of  cur vafcure  of  the 
sonic  line  Increases  in  che  region  rear  the  plane  of  sym¬ 
metry.  It  was  expected  that  the  sonic  line  would  terminate 
at  the  last  point  on  the  experimental  shock  curve.  As 
indicated  in  Fig.  7»  this  is  not  the  case.  Two  possible 
causes  of  this  difference  are  the  inriscid  flow  assumption 
and  error  in  the  numerical  technique. 

It  is  noted  that  the  first  family  characteristics  curve 
upward  as  region  III  is  traversed  in  the  inboard  direction. 
This  is  in  opposition  to  the  qualitative  sketch  of  the  cross- 
flow  plane  proposed  by  Bannink,  2fefebeling,  and  Reyn  and  sup¬ 
ports  the  views  of  Hays  and  Probstein  (Ref  9:532). 

The  computed  static  and  pitot  pressure  ratios,  P/P«, 
and  Pr/F  respectively,  at  £  -  0.266  are  given  in  Pig.  3. 
Ohis  figure  also  includes  exnerimental  data  for  P^/P_  taken 

P  Pcc 

from  Bannink ,  Kebbeling,  and  Reyn,  me  calculated  and 
experimental  pressure  ratios  agree  well  up  to  the  point  just 
prior  to  the  shock  wave  where  the  experimental  data  indi¬ 
cates  recompression  begins. 

Figure  9  is  a  plot  of  the  slope  of  the  streamline 
traces  in  the  crossflow  plane,  me  point  for  which  the 
streamline  slope  is  computed  is  at  the  center  of  each  trace. 
It  is  noted  that  the  streamline  traces  along  the  sonic  line 
are  nearly  parallel  to  the  soric  line.  However,  computa¬ 
tions  show  that  the  slope  of  the  sonic  line  is  slightly 
greater  chan  che  slope  of  the  streamline  traces. 
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Figure  8.  Pressure  Ratios  at  $=0,266  for  t*t0=2.96  a=l4.2 
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Iht-  i®t*i  Kach  nartb-g-r  distribution  st  {  =  0.256  is 
given  ia  rig.  13-  Uhe  kaoh  srziber  increases  uniformly 
from:  a  free  stream  value  of  2.56  to  a  raise  of  3.79  at  the 
crossflow  sonic  line. 

Stefce  Set  2 

The  characteristic  network  in  region  I  am  III  is 
©lotted  in  Fig.  ii  for  this  data  set.  Ibis  data  set  is  also 
one  of  the  sets  of  conditions  for  which  Babayev  presented  a 
solution.  she  solution  of  the  flow  field  in  region  I  as 
given  here  and  that  presented  by  Babayev  are  in  complete 
agree  ns  ns  •  However  ,  there  is  consioeraole  difference  be¬ 
tween  the  two  solutions  for  region  III.  The  crossflow  sonic 
line  given  by  Babayev  is  also  plotted  in  Fig.  11.  The 
difference  in  the  location  of  this  line  is  a  result  of  the 
planar  shock  assumption.  Experimental  data  is  not  available 
for  this  data  set  but,  in  an  analogous  study,  Bannink, 
Kebbeling,  and  Heyn  did  not  confirm  the  existence  of  a 
planar  shock  wave. 

The  computed  static  and  pitot  pressure  ratios  at 
i=  0.266  for  this  input  data  are  given  in  Fig.  12.  Figure 
13  is  a  plot  of  the  streamline  traces  and  Fig.  l4  gives  the 
Kach  number  distribution  at  £  =  0.266. 

The  program  was  run  with  this  set  of  input  data  and 
with  ten  Prandtl-Keyer  expansion  rays  in  region  I.  The 
most  s^ cni f'ieant  difference  was  that  the  terminal  position 
of  the  sonic  line  was  shifted  slightly  inboard  and  lower. 
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Figure  ll.  Characteristic  Network  for  a*  5»0°»  X»60.0 
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Characteristic  Network  for  a*l2»0 
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Data  Set  J 

The  characteristic  nefcatorh  for  this  data  set  Is  gives 
Ira  Pig.  15.  Figure  16  Is  a  plot  of  the  static  and  pitot 
pressure  ratios  at  «  =  0,266  and  Pigs.  17  and  IS  are  plots 
of  the  streamline  traces  am  the  Kach  number  distribution 
at  i  -  0.265,  respectively.  Sals  data  set  vss  used  bgr 
Jfeslerc  arho  assumed  that  no  shoes  existed  Ira  the  flow  field. 

Additional  Cases 

An  attempt  seas  made  to  obtain  solutions  for  K*  =  6.0, 
a=  7.0°  and  x=  60.0°  vlth  1GC  Pranifcl-Xeyer  expansion 
rays  and  for  =  2.96,  a  =  l4„2°  and  X  =  45.0°  (Data  set  1) 
with  ten  Pramtl-Xeyer  expansion  rays.  The  procedure  used 
to  determine  the  location  of  reflection  of  the  first 
expansion  ray  diverged  and  no  results  »fere  obtained. 
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V.  Conclusions  and  3eeoa2eBf3tIons 


Conclusions 

The  following  conclusions  are  nade  from  this  study. 

1.  It  is  possible  to  obtain  a  solution  for  the  flow 
field  in  the  pseudo-elliptic  region  by  use  of  the  method  of 
characteristics . 

2.  The  predicted  pressure  ratios  in  the  pseudo- 
elliptic  region  were  in  reasonable  agreement  with  the 
available  experimental  data. 

3.  Numerical  solutions  were  not  obtained  for  two  sets 
of  input  conditions  as  a  result  of  a  deficiency  in  the  pro¬ 
cedure  used  to  locate  the  reflection  of  the  first  expansion 
ray. 

4.  The  accuracy  of  the  terminal  position  of  the  cross- 
flow  sonic  line  was  improved  with  the  Increase  in  the  number 
of  expansion  rays. 

5.  The  slope  of  the  streamlines  along  the  sonic  line  is 
very  near  to  the  slope  of  the  sonic  line.  This  indicates 
that  most  of  the  fluid  in  the  pseudo-elliptic  region  exits 
this  region  through  the  curved  shock  wave.  Therefore,  a 
significant  portion  of  the  flow  field  governed  by  elliptic 
equations  is  rotational. 
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Reco-nerviatiors 

The  following  items  are  recommended  as  areas  for 
further  study. 

1.  An  attempt  should  be  made  to  improve  the  accuracy 
of  the  solution  by  modifying  the  procedure  used  to  locate 
the  sonic  points. 

2.  A  study  should  be  undertaken  to  improve  the 
technique  for  locating  the  reflection  of  the  first  expansion 
ray.  This  would  broaden  the  usefulness  of  the  procedure  and 
allow  for  further  characterization  of  the  flow  field. 

3.  An  experimental  investigation  of  the  flow  field 
similar  to  that  conducted  by  Bannink,  Nebbeling,  and  Reyn 
but  for  a  wide  range  of  free  stream  Mach  numbers,  angles  of 
attack,  and  sweep  angles  would  provide  additional  data  for 
further  verification  of  analytical  solutions. 

4.  A  study  of  the  elliptic  regions  should  be  conducted 
using  the  procedure  developed  in  this  study  to  determine  the 
boundary  conditions. 

5.  An  analytical  study  should  be  undertaken  to  include 
the  effects  of  viscosity  and  a  curved  wing. 


GAM  'AE/72-6 


Bibliography 


1.  Babayev,  D.  A.  "Numerical  Solution  of  the  Problem  of 
Flov;  Around  the  Upper  Surface  of  a  Triangular  Wing  By 
a  Supersonic  Stream.  "  U.S.S.R.  Comp.  Math .  Math. 
Physics,  2:296-309  (1962). 

2.  Bannlnk,  W.  J. ,  C.  Nebbeling,  and  J.  W.  Reyn.  Investi¬ 
gation  of  the  Flow  Field  on  the  Expansion  Side  of  a 
Delta  Wing  with  Supersonic  leading  Edges.  Report  VTH- 
128.  Delft,  The  Netherlands:  Technological  University 
Delft,  August  1965* 

3.  Beeman,  E.  R.  and  S.  A.  Powers.  A  Method  For  Determin¬ 
ing  the  Complete  Flow  Field  Around  Conical  wings  at 
Supersonlc/Hyperscnlc  Speeds .  AIAA  Paper  69--64£.  New 
York:  American  Institute  of  Aeronautics  and  Astronaut¬ 
ics,  June  1969. 

4.  Bulakh,  B.  M,  and  J.  W.  Reyn.  "Remarks  on  the  Problem 
of  the  Delta-Shaped  Wing."  Jh  Appl.  Math .  Hech. , 

21:  207-210  (1967). 

5.  Busemann,  A.  "Druke  auf  kegelformige  Spltzen  bei 

Bewegung  mlt  liberschallgeschwindigkei  t .  "  Z.  augew. 
Math.  Hech. ,  9:  496-498  (1929). 

6.  Chaing,  C.  W.  and  R.  D.  Wagner,  Jr.  Analysis  of  Super¬ 
sonic  Conical  Flows .  NASA  TN  D-58S4.  Washington: 
National  Aeronautics  and  Space  Administration,  1970. 

7.  Cross,  E.  J.  and  W.  L.  Hahkey.  Investigation  of  the 
Leeward  Side  of  a  Delta  V.’l  ng  at  Hypersonic ‘Speeds. 

AIAA  Paper  63^S’73.  New  York:  American  Institute  of 
Aeronautics  and  Astronautics,  June  1968, 

8.  Fowell,  L.  R.  "Exact  and  Approximate  Solutions  for  the 

Supersonic  Delta  Wing. "  J.  of  the  Aeronautical 
Sciences ,  23 :  709-720,  7 7"$  TAuguso  1 9 j'o ) . 

9.  Hayes,  W.  D.  and  R.  F.  Probstein,  Hypersonic  Flow 
Theory  -  Volume  I  (second  edition).  New  York: 

Acedemic  Press  Inc.,  1966. 

10.  Kitowski,  J.  V.  Unpublished  papers. 


42 


GAM/AE/72-6 


11.  Kutler,  P.  and  H.  Lonaz.  A  Systematic  Development  of 
the  Supersonic  Flow  Fields  Over  and  behind  :.-;inrs  and 
Wing- Body  Conf igu rations  Us i ng  a  Shock-Caoturing 
Flnlte-Dlff erence  Approach.  AIAA  Paper  71-99.  2few 
York:  American  Institute  of  Aeronautics  and  Astro- 

!  nautlcs;  January,  1971. 

12.  Liepmann,  H.  V/.  and  A.  Roshko.  Elements  of  Gasdy- 
namlcs.  New  York:  John  Wiley  and  Sons,  Inc.,  1957. 

13.  Maslen,  S.  H.  Supersonic  Conical  Plows .  NACA  TN 
2651 .  Washington:  National  Advisory  Committee  for 
Aeronautics,  March  1962. 


GSZ/2£/?2-6 


Appendix  A 


Development  of  Velocity  Sanation 


The  development  of  the  velocity  equations  as  given 
here  was  suggested  by  Kitcwskl  (Ref  10), 

The  total  velocity  vector  at  any  point  in  region  I 
can  be  separated  into  the  components  normal  and  parallel  to 
the  leading  edge  of  the  wing.  Resolving  the?  :  components 
into  components  along  the  XYZ  coordinates  and  dividing  by 
7*0  results  in 


XL 


n  t 

—  cos  (a. -6)  cos  x  +  —  sin  x 

Voo  1  Vo 


(A-l) 


V 


— •  cos(a1-6)  sin x +  —  cos  x 
V«,  1  Vo 


(A-2) 


w  =  —  sln(ai-6)  (A-3) 

V  1 

*00 

where  is  the  angle  between  the  surface  of  the  wing  and 
*he  component  of  the  free  stream  velocity  normal  to  the 
leading  edge  and  6  is  the  angle  through  which  the  normal 
component  of  the  free  stream  is  deflected;  and  6  are 

given  by  equations  (13)  and  (9)  respectively. 

The  component  of  velocity  parallel  to  the  leading  edge 
is  unaff ected  by  the  expansion  process  and  is  given  by 
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Vt  =  ^  cos  p  (A-4) 

where  p  is  the  angle  between  the  leading  edge  of  the  wing 
and  the  free  stream  velocity  vector  and  is  given  by 
equation  (14). 

The  component  of  velocity  normal  to  the  leading  edge  is 

vn  »  V  (A-5) 

Using  isentropic  flow  relationships  this  becomes 


V_  =  M  Coo 
n  n 


2 


Y-l 

1+ - M 

2 


2 

n  J 


and,  since  -  M*,  sin  P 

n oo 


(A-6) 


V  =  M 
n  n 


Y-l  2  Y-l 
1+ - - 


2  2 


Cco 


2 

n 


(A-7) 


Thus , 


n 


V„ 


Moo 


Y-l  2 

1  +  -  Moo  - 

2 


Y-l  2 


Mm  cos  P 


Y-l  2 
1  +  -  M 

2 


I 


(A-8 ) 


n 
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r  2 


or 


n 

Vo 


v^.1  I? 

Y^*  ‘io 


1  Y-l  o 

_  d  +  —  Iff)  - 


Y'1  2  Q 

— -  cos  @ 
Y+l 


2  1  »-Y-l  o 

“  (1  +  *T*  Hn) 


Y+l 


From  equation  (7),  is  given  by 


and  equation  (A-9)  becomes 


n  _ 


2  1  Y-l  2.  Y-l  2  k- 

- 5  (1  +  -  ifc) - cos  p 

Y+l  C  2  Y+l 


which  reduces  to 


n 


V„ 


(Y+l)  fJo  Y+l 


Y-l 

+  ( - )  (1  -  cos2^  ) 


1  + 


--sp 

Y-l  f  Y+l 


(A-9) 


(A-10) 


(A-ll) 


(A- 12) 
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Substituting  equations  (A-t)  and  (A-12)  Into  equations 
(A-l),  (A-2),  and  (A-3)  and  then  using  b  as  defined  by 
equation  (15)  results  in  equations  (10),  (ll)#  and  (12) 
respect „vely. 
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